c  ---------------------------------------------------------------------------
c  CFL3D is a structured-grid, cell-centered, upwind-biased, Reynolds-averaged
c  Navier-Stokes (RANS) code. It can be run in parallel on multiple grid zones
c  with point-matched, patched, overset, or embedded connectivities. Both
c  multigrid and mesh sequencing are available in time-accurate or
c  steady-state modes.
c
c  Copyright 2001 United States Government as represented by the Administrator
c  of the National Aeronautics and Space Administration. All Rights Reserved.
c 
c  The CFL3D platform is licensed under the Apache License, Version 2.0 
c  (the "License"); you may not use this file except in compliance with the 
c  License. You may obtain a copy of the License at 
c  http://www.apache.org/licenses/LICENSE-2.0. 
c 
c  Unless required by applicable law or agreed to in writing, software 
c  distributed under the License is distributed on an "AS IS" BASIS, WITHOUT 
c  WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the 
c  License for the specific language governing permissions and limitations 
c  under the License.
c  ---------------------------------------------------------------------------
c
      subroutine getdhdr(datahdr,ibctyp,ndata)
c
c     $Id$
c
c***********************************************************************
c     Purpose:  To set character data for main output file headers 
c     when the 2000 series bc's are used 
c***********************************************************************
c
#   ifdef CMPLX
      implicit complex(a-h,o-z)
#   endif
c
      character*10 datahdr(11)
c
      do 1 m=1,10
    1 datahdr(m)  = '          '
c
      if (ibctyp.eq.2002) then
         datahdr(1)  = '    p/pinf'
      end if
c
      if (ibctyp.eq.2003) then
         datahdr(1)  = '      mach'
         datahdr(2)  = '   pt/pinf'
         datahdr(3)  = '   tt/tinf'
         datahdr(4)  = '     alpha'
         datahdr(5)  = '      beta'
         if (abs(ndata).gt.5) then
            datahdr(6)  = '     turb1'
         end if
         if (abs(ndata).gt.6) then
            datahdr(7)  = '     turb2'
         end if
      end if
c
      if (ibctyp.eq.2009 .or. ibctyp.eq.2010) then
         datahdr(1)  = '   pt/pinf'
         datahdr(2)  = '   tt/tinf'
         datahdr(3)  = '     alpha'
         datahdr(4)  = '      beta'
         if (abs(ndata).gt.4) then
            datahdr(5)  = '     turb1'
         end if
         if (abs(ndata).gt.5) then
            datahdr(6)  = '     turb2'
         end if
      end if
c
      if (ibctyp.eq.2019) then
         datahdr(1)  = '  pt/ptinf'
         datahdr(2)  = '  tt/ttinf'
         if (abs(ndata).gt.2) then
            datahdr(3)  = '     turb1'
         end if
         if (abs(ndata).gt.3) then
            datahdr(4)  = '     turb2'
         end if
      end if
c
      if (abs(ibctyp).eq.2004) then
         datahdr(1)  = '   tw/tinf'
         datahdr(2)  = '        cq'
      end if
c
      if (abs(ibctyp).eq.2014) then
         datahdr(1)  = '   tw/tinf'
         datahdr(2)  = '        cq'
         datahdr(3)  = '     index'
      end if
c
      if (abs(ibctyp).eq.2024) then
         datahdr(1)  = '   tw/tinf'
         datahdr(2)  = '        cq'
         datahdr(3)  = ' intermtbc'
      end if
c
      if (ibctyp.eq.2005) then
         datahdr(1)  = '      nblp'
         datahdr(2)  = '     dthtx'
         datahdr(3)  = '     dthty'
         datahdr(4)  = '     dthtz'
      end if
c
      if (ibctyp.eq.2006) then
         datahdr(1)  = '      nblc'
         datahdr(2)  = '    p/pinf'
         datahdr(3)  = ' integ dir'
         datahdr(4)  = ' axial dir'
      end if
c
      if (ibctyp.eq.2007) then
         datahdr(1)  = '       rho'
         datahdr(2)  = '         u'
         datahdr(3)  = '         v'
         datahdr(4)  = '         w'
         datahdr(5)  = '         p'
         if (abs(ndata).gt.5) then
            datahdr(6)  = '     turb1'
         end if
         if (abs(ndata).gt.6) then
            datahdr(7)  = '     turb2'
         end if
      end if
c
      if (ibctyp.eq.2102) then
         datahdr(1) = '    p/pinf'
         datahdr(2) = '   dp/pinf'
         datahdr(3) = '    rfreqp'
         datahdr(4) = '      lref'
      end if
c
      if (ibctyp.eq.2103) then
         datahdr(1) = '    p/pinf'
         datahdr(2) = '   dp/pinf'
         datahdr(3) = '    rfreqp'
         datahdr(4) = '      lref'
         datahdr(5) = '    phioff'
      end if
c
      if (ibctyp.eq.2008 .or. ibctyp.eq.2038) then
         datahdr(1)  = '       rho'
         datahdr(2)  = '         u'
         datahdr(3)  = '         v'
         datahdr(4)  = '         w'
         if (abs(ndata).gt.4) then
            datahdr(5)  = '     turb1'
         end if
         if (abs(ndata).gt.5) then
            datahdr(6)  = '     turb2'
         end if
      end if
c
      if (abs(ibctyp).eq.2016) then
         datahdr(1)  = '   tw/tinf'
         datahdr(2)  = '        cq'
         datahdr(3)  = '       cqu'
         datahdr(4)  = '     sjetx'
         datahdr(5)  = '     sjety'
         datahdr(6)  = '     sjetz'
         datahdr(7)  = '     rfreq'
      end if
c
      if (abs(ibctyp).eq.2026) then
         datahdr(1)  = '      vmag'
         datahdr(2)  = '     rfreq'
         datahdr(3)  = '  sideangj'
         datahdr(4)  = '       sxa'
         datahdr(5)  = '       sya'
         datahdr(6)  = '       sza'
         datahdr(7)  = '       sxb'
         datahdr(8)  = '       syb'
         datahdr(9)  = '       szb'
         if (abs(ndata).gt.9) then
            datahdr(10)  = '     turb1'
         end if
         if (abs(ndata).gt.10) then
            datahdr(11)  = '     turb2'
         end if
      end if
c
      if (ibctyp.eq.2018) then
         datahdr(1)  = '         T'
         datahdr(2)  = '      rhou'
         datahdr(3)  = '      rhov'
         datahdr(4)  = '      rhow'
         if (abs(ndata).gt.4) then
            datahdr(5)  = '     turb1'
         end if
         if (abs(ndata).gt.5) then
            datahdr(6)  = '     turb2'
         end if
      end if
c
      if (ibctyp.eq.2028) then
         datahdr(1)  = '      freq'
         datahdr(2)  = '   rhoumax'
         datahdr(3)  = '   rhovmax'
         datahdr(4)  = '   rhowmax'
         if (abs(ndata).gt.4) then
            datahdr(5)  = '     turb1'
         end if
         if (abs(ndata).gt.5) then
            datahdr(6)  = '     turb2'
         end if
      end if
c
      return
      end
